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ABSTRACT 

Using a full spectral scan of an active region from the Extreme-Ultraviolet 
Imaging Spectrometer (EIS) we have obtained Emission Measure EM(T) distri- 
butions in two different moss regions within the same active region. We have 
compared these with theoretical transition region EMs derived for three limiting 
cases, namely static eq uilibrium, strong condensation and strong evaporation from 



Klimchuk et al.l ( 120081 ) . The EM distributions in both the moss regions are strik- 
ingly similar and show a monotonically increasing trend from log T[K] = 5.15-6.3. 
Using photospheric abundances we obtain a consistent EM distribution for all 
ions. Comparing the observed and theoretical EM distributions, we find that 
the observed EM distribution is best explained by the strong condensation case 
(EM con ), suggesting that a downward enthalpy flux plays an important and pos- 
sibly dominant role in powering the transition region moss emission. The down- 
flows could be due to unresolved coronal plasma that is cooling and draining 
after having been impulsively heated. This supports the idea that the hot loops 
(with temperatures of 3-5 MK) seen in the core of active regions are heated by 
nanoflares. 



Subject headings: Sun: corona — Sun: atmosphere — Sun: transition region - 
Sun: UV radiation 
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Introduction 



Solar coronal heating has been one of the most intensely studied problems in solar 
physics over the past 60 years. Despite ma jor advances in t he observational and theoretical 
tools, the solution remains elusive (see e.g.. lKlimchuk II2006I ). Active regions provide a target 
of opportunity to study the dominant heating mechanisms for the corona. The basic building 
blocks of active regions are coronal loops. Because of the lack of cross-field transport of mass 
and thermal energy in the corona, the problem of coronal heating can effectively be reduced 
to a clear understanding of the heating mechanism in a single isolated coronal loop or loop 
strand. 

The large-scale 1 MK loops and fan-like loops (which appear at the periphery of active re- 
gions) are seen c learly in observation s recorded using the Transition Region and Coronal Ex- 
plorer (TRACE; lHandy et al. Ill999l ) 171 A passband. Their properties appear to be consis- 

tent with impulsive heating models (see e.g., Tripathi et al. 112009 ; Warren. Winebarger. and Mariska 
20031 ; Klimchuk 112006k 12009k and references cited therein). However, the 3-5 MK loops in the 
core of the active regions are more diffuse and more difficult to isolate even with present day 
observatories. Hence it is more difficult to measure the physical plasma parameters along 
these hot core loops and to compar e those properties with simulated results (although see 
Lopez Fuentes and K limc huk 1120101). The "moss" regions are located at footpoi nts of hot 



(3-5 MK) loops (Martens et al. 



2000; Antiochos et al. 



map the photospheric magnetic flux very closely ( Tripathi et al. 20081 ). The moss regions 



20031: iTripathi et al.llioioh and they 



being the footpoints of the hot loops, provide a unique opportunity to study the dominant 
heating mechanism at work in active region core loops. 

A number of publications involving observational results, analytical and ID numerical 
simulations have proposed that steady heating is at play in the hot loops associated with 
moss regions. The main argument is that minimal variability is observed in brightness of the 



moss emission or in the Doppler shifts and line widths of spectral lines. lAntiochos et al. 



( 120031 ) pointed out that impulsive heating would produce transient loops in the cores of 
active regions that would be seen by TRACE as the plasma cools through 1 MK. Such 
loops are present but not ubiquitous. The authors also no ted that the brightne s s of m oss is 
quite steady when averaged over several TRACE pixels. iBrooks and Warren I ( 120091 ) later 
showed that Doppler shifts and line widths observed by t he Extreme ultraviole t Imaging 
Spectrometer (EIS) on Hinode are also quite steady. Finally, ITripathi et al.l ( 120101 ) measured 
the density and thermal properties of moss observed by EIS and found that they do not 
change significantly over periods ranging from an hour to 5 days. 



These observations provide compelling evidence that heating in the cores of active re- 
gions must be quasi-steady if the cross- field scale of the heating is comparable to the 2 arcsec 
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(1.5 Mm) resolution of EIS or larger. That is, we can conclude that the heating is steady, but 
only if the plasma does not have unresolved sub- structure. There is good reason to believe 
that such sub-structure may exist, however (e.g.. iKlimchuk 1 120061 ) . Individual thin strands 



could be highly time variable, and yet an unresolved bundle would appear steady as long as 
the strands are not in phase. This is true for both the coronal portions of the strands and 
their moss footpoints. 

We therefore seek a diagnostic of heating that does not require that the fundamen- 
tal structures be resolved. Doppler shifts are one possibility. Steady heating tends to 
produce no Doppler shif t because the loop evo lves to a static equilibrium (although see 

20101 ). Impulsive heating would produce a net 



Mariska and Boris 1983; Klimchuk et al. 



red shift on the other hand. In an unresolved bundle of strands, some strands would have 
upflows due to evaporation and others would have downflows due to condensation (cooling 
and draining). Simulations show that the upflows are faster, fainter, and shorter lived than 
the downflows (IPatsourakos and Klimchuk 1 120061 ) . Hence, the composite line profile would 
be red shifted by several kilometers per second, and very hot lines would also have a faint 
blue wing enhancement. The red shifts are pre dicted to decrease with temper ature and 
may even become blue shifts in the hottest lines ( jBradshaw and Klimchuk Il2010l in prepa- 
ration). The observational picture concerning Doppler shifts in moss regions is unclear. 
Some published measurements indicate minimal red-shifts consistent with steady heating 



Brooks and Warren 1120091 ). w hile other measurements indicate larger red-shifts con- 



le.g., 

sistent with impulsive heating (e.g ., 



Doschek et al.lll97a iBrueckner Ill98ll : Klimchuk I Il987l : 



Doschek et al. Il2008t iDel Zanna 1120081 ). This is an important area of future work. 



Another diagnostic of coronal heating that works when there is unresolved sub-structure 
is the emission measure distribution, EM(T), or the closely related differential emission mea- 
sure distribution, DEM(T). Since the plasma is roughly isothermal along the coronal portion 
of a strand, EM(T) in the corona is determined primarily by the distribution of strand tem- 
perature s. An isothermal corona with all str ands having the same temperature implies steady 
heating ( Warren 1999; Eiandi et al~ 2002), but some coronal EM(T) suppo r t a nanoflare 



interp r etation (e.g., 



20091 ). IWarren et al. 



atsourakos and Klimchuk 1 120091 ; lO'Dwyer et al. I I201CH : iReale et al. 



( 201o[ ) suggested that the lack of a strong correlation between hot and 
warm emission in the core of an active region is evidence that the heating must be steady. 
We disagree that this is the only interpretation. The relative intensities of coronal emission 
formed at, say, logT = 6.0 and 6.5 can vary greatly depending on the magnitude of the 
nanoflares (e.g.. IKlimchuk et al.ll2008l Figs. 4 and 7, dashed curves). 



Moss is the transition region footpoints of hot strands, an d the EM(T) there is deter- 
mined largely by the variation of temperature along the strands. IKlimchuk et al.l (120081 ) have 
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shown that the form of EM(T) is different under three limiting conditions: static equilibrium, 
strong condensation and strong evaporation. We make use of these differences in the work 
reported here. We measure the EM distribution of two moss regions observed by EIS and 
compare them to the theoretical predictions. In this way we distinguish between steady and 
impulsive heating. 

The rest of the paper is structured as follows. In section [2] we p resent our obser v ations . 
In section [3] we provide the emission measure expressions following iKlimchuk et al.l ( 120081 ). 
Section H] provides results and we summarise and discuss our results in section 



2. Observations 



In this study, we have analysed observations recorded by the EIS aboard Hinode. The 
EIS instrument observes the spectral ranges 170-210 and 250-290 A, providing observations 
in a broad ran g e of t emperatures (logT[K] = 4.7-7.3). For detailed information on EIS see 
Culhane et al.l ( 120071 ). For most of the EIS studies only a limited number of spectral lines 
were chosen due to the limitations on telemetry. However, occasionally full spectral scans of 
specific regions have been telemetered. Here we use one such raster of an active region AR 
10961 recorded on July 01, 2007 at 03:18 UT. The active region first appeared on the east 
limb on June 26th, 2007. The data analysed were obtained from a single raster that scanned 
a field of view (FOV) of 128" by 128" with the 1" slit, moving from west to east with an 
exposure time of 25 sec. The over-plotted box on the TRACE 171 A image displayed in the 
left panel of Fig. [1] marks the region which was rastered by EIS. An EIS image obtained in 
Fe XII A195.12 A is shown in the right panel of Fig. [TJ 

The full spectral scan allowed us to choose spectral lines formed over a broad range of 
temperature, which is essential in order to obtain a good EM distribution of the plasma. 
This is the first time a full spectral scan is being used to derive an EM distribution in active 
region moss. In our study we have used relatively clean spectral lines (see Table ED) with tem- 
peratures of formation logT[K ] = 5.1-6.5 based on EIS spectral atlases (ILandi and Young 



20091 : lYoung and Landi lbo09h . The standard EIS software provided in SolarSoft were used 
to process the data and eis-autoJit, also provided in SolarSoft, was used for Gaussian fitting 
the spectral lines. To de-b lend some of th e blend ed lines labelled 'b' in Table H] we have used 



the same pr ocedure as in iTripathi et al.l (120101 ) based on the suggestions of 
( l2007Ll2009f ). 



Young et al. 



We have derived EM following the approach of iPottasch I (119631 ) whereby individual 
emission lines yield estimates of the emission measure at the lines' temperatures of formation. 
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Table 1: Spectral lines used to study the emission measure distribution in moss regions. 



Ion 


Wavelength 


log T 




[A] 


[K] 


IV 


279.6 


5.15 


IV 


279.9 


5.15 


V 


248.5 


5.35 


VI 


183.9 


5.45 


VI 


184.1 


5.45 


Mg V 


276.6 


5.45 


Si VI 


246.0 


5.60 


Fe VIII 


185.2 


5.60 


Mg VI 


269.0 


5.65 


Mg VII 6 


278.4 


5.80 


Si VII 


275.3 


5.80 


Fe IX 


171.1 


5.85 


Fe X 


184.5 


6.05 


Si IX 


258.1 


6.05 


SiX 


261.0 


6.15 


Fe XI 


180.4 


6.15 


Fe XI 


188.3 


6.15 


Fe XII 


192.4 


6.20 


Fe XII 


193.5 


6.20 


Fe XII 6 


195.1 


6.20 


Fe XIII 


202.0 


6.25 


Fe XIV 6 


274.2 


6.30 


Fe XV 


284.2 


6.35 


Fe XVI 


263.0 


6.45 


Ca XIV 


193.8 


6.55 
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An EM distribution is built up by considering lines formed over a wide range of temperatures. 
The method requires the contribution function to be approximated by a function such that 
it is defined to be a constant over the temperature range log T max — 0.15 to log T max + 0.15 
where T max is the temperature where the contribution funct ion has its maximum. A detailed 
description of this method is given in iTripathi et al.l (120101 ) . 



3. Theoretical Transition Region Emission Measure 



Klimchuk et al.l (120081 ) derived the differential emission measure (DEM) distribution of 
the transition region of a strand in the limiting cases of strong evaporation (Equation A3), 
strong condensation (Equation A7), and static equilibrium (Equation A12). The transition 
region is taken here to be the region of steep gradient at the strand footpoints. Its tem- 
perature ranges from a chromospheric value up to approximately one-half of the maximum 
temperature in the strand. Thus, it can reach several MK for strands which are very hot. 
In the case of static equilibrium, a downward heat flux from the corona powers the radiative 
losses from the transition region. These losses are ignorable in the case of strong evaporation, 
where a downward heat flux drives an evaporative upflow. The energy flux of the upflow is 
dominated by enthalpy if the flow is subsonic. For the final case of strong condensation, the 
enthalpy of a downflow powers the t ransi tion r egion radiation, with the heat fl ux being rela- 
tively unimportant. See iBradshawl (120081 ) and lBradshaw and Cargill for detailed 
discussions of the importance of downflows in heating the transition region and cooling the 
corona. 



The expressions given in lKlimchuk et al.l ( 120081 ) are for the differential emission measure. 
The conversion to emission measure (EM) requires that they be multiplied by temperature 
and by Zn(10), giving 



EM se ps Zn(10) 



1/2 p ^3/4 

k A(T)V2 



... static equilibrium 



5 k J T 

EM con ps — Zn(10) — a(T) — " s ^ ron ^ condensation 



(2) 



EM ev ps 



20 ¥ Jo 



strong evaporation 



(3) 



where Kq— 1.0x10 6 in cgs units, k is Boltzmann's constant, T is temperature. Jo is mass 
flux defined as Jo — nv , where n and v are electron number density and plasma flow speed 
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respectively. P is average pressure along a strand. A(T) is the optically thin radiative 
loss function. Fig . [2] displays the radiative losses as a function of temperature calculated 
with CHIANTI ( IDere et al.l Il997l ) by free-free, radiative r ecombination and by line r adia- 
tion for photospheric abundances (dashed-dotted line) of iGrevesse and Sauval and 
coronal abundances (solid line) of Feldman"! (1992 ). For t his calculation, we have used the 
CHIANTI v6.0.1 ionization equilibrium (IDere et al.l 120091 ) . The above equations give the 
emission measure curves for individual loop strands. While computing the theoretical EMs 
to compare with the observed EMs, we have considered the average pressure (P) and mass 
flux Jo as arbitrary constants. We have restricted ourselves to the temperature dependence 
of EM e „, EM con and EM se and compared the shapes of the theoretical EM curves with those 
derived observationally as explained in section [2j The assumption of constant pressure and 
constant mass flux in a given strand is justified so long as the different parts of the solar 
atmosphere are connected along the magnetic field. However, this assumption would break 
down if plasma in differe nt layers of the atmo s phere were physically in different structures 
as was suggested by e.g., lLandi and Feldman I ( 120081 ) and references therein. 



Results 



Figure E] displays EM curves for two different moss regions labelled "A" and "B" in 



the right pane 
(top pane l ) of 



of Fig. [Tj The EM curves are created using the photospheric abundances 



Grevesse and Sauvall (119981 ) and the coronal abundances (bottom panel) of 



Feldman I (119921 ). The error bars on the EMs are computed by considering a 30% error in 
the observed intensities for all the spectral lines, which is a rough estimate of the cumulative 
errors due to radiometric calibration, line fitting and atomic calculations. 

The EM curves for the two different regions displayed in Fig. [3] are very similar to each 
other when the same abundances are used. However, for a given region the EM curves 
are significantly different for the two different abundances, mainly at lower temperatures. 
The difference in the trend of the EM curves is basically due to the three data points 
on the extreme left, corresponding to three Oxygen lines namely O IV, O V and O VI. 
The EMs for Oxygen lines do not appear to be consistent with the EMs for other lines 
when coronal abundances are used. However, this discrepancy disappears when we use 
photospheric abundances. From the figure we note that the EMs of the Oxygen lines do 
not change when coronal or photospheric abundances are used. However, for other lines, the 
EMs increased by a factor of four for photospheric abundances, which makes the EMs for 
other lines consistent with the EMs for the Oxygen lines. Therefore, it seems appropriate to 
conclude that the composition in these moss regions is photospheric. However, we note that 
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Oxygen is the only element in our dataset with a high first ionization potential (FIP). The 
high-FIP O VI lines agree with the low-FIP Mg V line at logT = 5.4 only with photospheric 
abundances. To verify the conclusion that moss regions must have photospheric abundances, 
it is important to examine additional pairs of high-FIP and low-FIP lines formed at the same 
temperature. 

The o bservational picture of abund ances in the solar upper atmosphere is not well estab- 



lished (see iFeldman and Widing 1120031 for a recent review). Using observations recorded by 



Skylab S082A, it was found that in newly emerging active regions the composition was nearly 
photospheric. However, within a few hours after their emerg ence the composition starts to 
modify and it approaches a typical coronal composition (see ISheeley I Il995l ; IWiding 1 119971 ; 
Widing and Feldman Il200l[ ). Using data from the Corona l Diagnostic Spectrometer (CDS) 
aboard the Solar and Heliospheric Observatory (SoHO), lYoung and Mason I (119971 ) found 
photospheric composition in an emerging flux region. The active region studied here had 
emerged about 6 days prior to the date when the data analysed here were recorded. Consid- 
ering the results obtained using the Skylab observations, our analysis should suggest coronal 
abundances. It is worth mentioning, however, that Skylab S082A was a slitless spectrom- 
eter whi ch produced overl apping spectroheliorams providing partial information for active 



regions. iDel Zanna I (120031 ) suggested that the previous results based on Skylab observations 
overestimated the FIP bias by a large factor for active regions, if the plasma is isothermal 
rather than multi-thermal. Furthermore, since the moss regions are presumably confined to 
a relatively low altitude in the solar atmosphere, the FIP effect might not be significant in 
these regions. Therefore, photospheric abundances could prevail in moss regions. 

The EM curves for the moss regio ns show that the pea k of the emission comes from 
Fe XIII (logT = 6.25) as suggested bv iTripathi et all feoioh . In addition, the EM(T) ob- 
tained using photospheric abundances shows a monotonic increase at lower temperatures 
from logT = 5.15. However, the EM(T) obtained using coronal abundances shows a min- 
imum in the emission measure at ~ logT = 5.5 similar to the EM curve obtained us- 
i ng the DEM distribution derived from SERTS EUV spectrum using the same abundances 
dBrosius et al~lll996[ lioooh . 



We have compared the observed EM distributions obtained using coronal and photo- 
spheric abundances with those predicted theoretically for the three limiting cases of static 
equilibrium (EM se , Eq. QP, strong condensation (EM con , Eq. [2]) and strong evaporation (EM e „, 
Eq. [3]) for both the regions 'A' and 'B' marked in the right panel of Fig. [TJ Fig. H] displays 
the observational EMs obtained using photospheric abundances over-plotted with the the- 
oretical EMs for the three limiting cases, both for region 'A' and region 'B' in the left and 
right columns respectively. 
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The theoretical curve that provides the best overall match with the observed points is 
EM con i.e., the strong condensation case (see plots in the middle row in Fig. H|). The EM se 
for static equilibrium matches quite well with observed EM for temperatures below 1 MK 
but not at higher temperatures. However, that for strong condensation (EM con ) matches 
quite well for logT = 5.2 - 6.25. Above logT = 6.3, the EM con suggests an increasing 
trend, whereas the observational EM starts to decrease. We can understand this apparent 
discrepancy as follows. 

Suppose that unresolved strands are heating impulsively to very high temperatures by 
nanoflares and subsequently cooling. As stated earlier, the transition region in a given strand 
extends to roughly one-half of the maximum strand temperature. Recently heated strands 
will therefore contribute to moss emission over a wide range of temperatures, while strands 
that have cooled will only contribute emission at lower temperatures. We can imagine a 
collection of curves like those in the middle panels of Fig. H] except truncated at different 
temperatures above logT = 6.3. The sum of those curves would have a shape similar to 
the observed EM(T), rolling over at logT = 6.3. This basic effect can be seen in the 



transition region of model 4 in iKlimchuk et al.l (120081 Fig. 7, dot-dashed curve). The figure 



shows DEM(T) (= EM(T)/T) and is therefore flat rather than positively sloped at the lower 
temperatures. The temperature where EM(T) is expected to roll over will depend upon the 
number distribution and brightness of strands in different stages of cooling. Strands which 
have cooled below 1 MK should be relatively faint due to mass draining and therefore will 
have a reduced contribution to the composite EM(T). 

The EM ev curve for strong evaporation does not show any resemblance to the observed 
EM curve. This is not surprising, since, in the nanoflare scenario, evaporative upflows 
are much fainter and much shorter lived than the subsequent condensation downflows and 
therefore contribute minimally to the net emission. 

When using coronal abundances, we did not find a reasonable correspondence between 
the observed EM(T) and any of the theoretical EM(T). However, this conclusion is based on 
the disagreement in the emission measure at logT < 5.45, and therefore depends entirely on 
EM of Oxygen lines. The EM con and EM se curves fit the data equally well in this questionable 
range logT < 5.45. However, we note that at logT > 5.45 EM con provides a better fit to 
the data than EM se regardless of which abundance set is assumed. 
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5. Summary and discussion 



We have obtained EM distributions in active region moss using full spectral scans 
recorded with EIS for two different regions. The EM distributions for the two regions 
show very good agreement with each other. The EM obtained using coronal abundances 
does not give consistent values for different lines. However, photospheric abundances seem 
to provide consistent EMs. The EM distributions show a monotonically increasing trend 
from logT = 5.15 - 6.3. Of the three theoretical limiting cases we have considered, strong 
condensation (EM con ) best reproduces the observations. Static equilibrium (EM con ) does a 
good job at lower temperatures, but it fails significantly above logT = 6.0. Strong evap- 
oration (EM con ) is clearly inconsistent with the observations. We conclude from this that 



Bradshawl (120081) on the basis of observed red-shifts (see e.g. .iDoschek et al. 



reached by 


Bradshaw 


(2008) 


Brueckner 


1981 




Klimchuk 


2008 


; iTripathi et al. 


2009). 



1987; Doschek et al. 2008; Del Zanna 2008; Warren et al 



1976 



As we have discussed, one obvious source of downflows is the cooling and draining of 
plasma that has been impulsively heated by nanoflares. As long as the heated strands have 
a sub-resolution diameter, this picture is consistent with observations of apparently steady 
intensities, Doppler shifts, and line widths. Another possible source of downflows is material 
from Type II spicules. This newly discovered phenomena appears to rep resent jets of chro 



mosp heric plasma that is heated to coronal temperatures as it is ejected (IDe Pontieu et al 



20071 ). Presumably the material subsequently cools and falls. Much more observational and 
modeling work is necessary before any firm conclusions can be drawn. 

Although our results support the picture of nanoflares and possibly Type II spicules in 
the cores of active regions, we cannot rule out steady heating entirely. If magnetic flux strands 
constrict appreciabl y at their base, then EM se (T) will have a different shape from that which 
we have shown (e.g. IWarren et al. 1120101 ) . The theoretical curve could in principle be brought 
in line with observations. However, because the transition region is so thin, the constriction 
must be very dramatic and the transition region must fall right where the constriction occurs. 
This would seem rather surprising given that the height of the transition region depends on 
the coronal pressure and moves up and down in response to changes in the heating rate, either 
slow or rapid. It has been found that the area of qui et Sun network boundari es increases 
by a factor of about 1.6 between logT = 5.4 and 6.0 (IPatsourakos et al. Ill999l ). How this 
relates to moss regions discussed here is unclear. 

The result that photospheric abundances apply in regions of downflow may hold impor- 
tant clues about the origin of the FIP effect, the fractionation of elements based on their 
first ionization potential. If fractionation (coronal abundances) occurs in regions of steady 
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heating, where conditions are static, but not in regions of impulsive heating, where there are 
continual upflows and downflows, then it is suggested that the fractionation process involves 
some sort of long-term settling. Further observational and modeling work is needed also in 
this regard. 
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Fig. 1. — Left panel: TRACE 171 A image showing the complete active region. The over- 
plotted box shows the EIS field of view. Right panel: An image obtained in Fe XII 195.12 A. 
The two regions labelled 'A' and 'B' were selected for detailed study. 




Fig. 2. — Radiative energy losses using the photospheric abunda nces (dash e d-dot ted line) of 
Grevesse and Sauvall (119981 ). coro nal abundances (s olid line) of iFeldman I (119921 ) and CHI- 
ANTI v6.0 ionization equilibrium iDere et al.l (120091 ) . 
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Fig. 3. — Observed EM(T) for two different moss regions A and B shown in the right p anel 
of Fig. [1] using the photospheric abundanc es (top pa n el) of iGrevesse and Sauvall (119981 ) and 
the coronal abundances (bottom panel) of iFeldman I (11992h . 
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Fig. 4. — Theoretical (solid lines with over-plotted diamonds) EM curve for three limiting 
ca ses using Eqns .1, 2 fc 3 and observed (asterisks) EM curves using photospheric abundances 
of iGrevesse and Sauvall (119981 ). 



